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We derive an effective quantum Josephson array model for a weakly interacting one-dimensional 
condensate that is fragmented into weakly coupled puddles by a disorder potential. The distribution 
of coupling constants, obtained from first principles, indicate that weakly interacting bosons in a 
disorder potential undergo a superfluid insulator transition controlled by a strong randomness fixed 
point [Phys. Rev. Lett. 93, 150402 (2004)]. We compute renormalization group flows for concrete 
realizations of the disorder potential to facilitate finite size scaling of experimental results and allow 
comparison to the behavior dictated by the strong randomness fixed point. The phase diagram of 
the system is obtained with corrections to mean-field results. 



I. INTRODUCTION 
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A number of recent experiments have investigated the 
properties of Bose condensates in disordered traps ^ and 
thereby revived the theoretical interest^^ in a fundamen- 
tal unsolved problem: what is the fate of Anderson local- 
ization in the presence of interactions and strong quan- 
tum correlations? The problem is particularly intrigu- 
ing in one dimension where disorder alone or interactions 
alone would have a profound effect on the physics. 

One way to tackle the problem is to start from the 
harmonic fluid (Luttinger liquid) description of a uni- 
form interacting Bose gas^, then add disorder to it as a 
small perturbation. Using this approach Giamarchi and 
Schultz^ predicted a superfluid-insulator quantum phase 
transition that occurs at a universal value of the Lut- 
tinger parameter, or correlation decay exponent. This 
approach is justified a priori if the chemical potential 
set by the interactions is much larger than the disorder 
strength. However, most experiments are in the oppo- 
site limit of weak interactions, where the disorder acts to 
fragment the condensate invalidating the Luttinger liq- 
uid description. The non-interacting ground state is also 
a bad starting point. This state, in which all particles oc- 
cupy the lowest single-particle localized state, is unstable 
to adding even the weakest interaction. The absence of 
a simple basis from which to formulate a perturbation 
expansion makes the limit of weakly interacting bosons 
inherently strongly correlated. 

In this paper we derive from first principles a low en- 
ergy effective model of the disordered quantum gases in 
the form of a random Josephson junction array. Such a 
model was assumed in previous work, coauthored by one 



of 



^,4,10,11 



as a starting point for a real-space renor- 



malization group (RSRG) analysis^ This analysis 
predicts a superfluid insulator transition controlled by 
a strong randomness fixed point, distinct from the tran- 
sition described by Giamarchi and Schulz^. However the 
connection between the quantum gas in the disorder po- 
tential and the effective Josephson array model has not 
been established. 

Starting from the microscopic random potential for the 
atoms, we show how the condensate fragments into mezo- 
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FIG. 1: Schematic sketch of the disorder potential and conden- 
sate fragmentation, (a) Regime of smooth disorder - typical 
puddles are much larger than the healing length, (b) A rough 
disorder is effectively smoothed by the healing length or single 
particle localization length. In both cases the distribution of 
weak links is determined by atypical long barriers. 



scopic puddles. We derive the distribution of Josephson 
couplings between the puddles and of the charging ener- 
gies within them. This is done for two different regimes 
of the disorder potential: rough and smooth potential as 
compared to the healing length and to the single particle 
localization length of the condensate (see Fig. 1). Re- 
markably in both cases the distributions calculated from 
first principles, are precisely in the form of the stable 
solutions of the RG equations found in Ref. 4. In par- 
ticular the distribution of Josephson links P{J) behaves 
as a power law of J at small J. Hence given a set of mi- 
croscopic parameters we can immediately compute the 
phase diagram and make direct predictions for finite size 
scaling of observables. 
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II. THE MODEL AND MAPPING TO 
JOSEPHSON ARRAY MODEL 

Our starting point for the theoretical analysis is the 
continuum boson hamiltonian 



H 



dxip"^ [ Vix) — /i I + m/j^t/j^ipilj (1) 

\ 2m J 



Here u is the effective contact interaction and V{x) is 
a random potential assumed to be gaussian and charac- 
terized by the auto-correlation function {V{x)V{x^)) ^ 

For a given chemical potential, the potential landscape 
is filled with particles up to the chemical potential, form- 
ing local superfluid puddles (see Fig. 1(a)). Each puddle 
is assumed to be characterized by a single phase (pi and 
occupation number Ni which are non-commutating op- 
erators. We would like to map the physics of this system 
to an effective random Josephson junction array model 



Here Ui is the inverse capacitance, or charging energy of 
a of puddle i and Ji is the Josephson coupling between 
puddles i and i -\- 1. 

The Josephson couplings can be obtained from the 
mean field Gross-Pitaevskii (GP) solution of (1). The 
formula for the superfluid stiffness of a one-dimensional 
condensate in terms of the GP wavefunction is (see 
appendix A and Ref. 14) 
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This integral is dominated by deep minima of that oc- 
cur between neighboring puddles. We can compute the 
stiffness as a sum over those saddle point contributions 
Paif ~ l/^spX]^r^' where the Ji are interpreted as 
the Josephson couplings between puddles and Nsp is the 
number of deep minima. More generally, the integral (3) 
can be subdivided into sub regions of order of the corre- 
lation of |V^(x)p to obtain a distribution of local phase 
stiffness. The same analysis can be used to obtain the ef- 
fective Josephson array from a density profile measured 
in experiment. We emphasize that this is a method to 
obtain the local phase stiffness Ji between puddles, which 
enter the effective Hamiltonian (2). The true thermody- 
namic stiffness can then be calculated only within the 
quantum Hamiltonian (2). In general it will be renor- 
malized downward compared to the mean field stiffness 
(3) because of phase slips induced by the charging terms 

We shall compute ip{x) using a simple approximation 
to the GP equation, which allows us to obtain analytic 
results for the distributions of coupling constants. The 
results will be checked against the distributions obtained 
by exact numerical calculation of the GP ground state. 



We distinguish two regimes according to the ratio be- 
tween the length scale a of disorder potential fluctuations 
to the natural correlation length on which the conden- 
sate amplitude can adopt to the external changes. 
The latter is determined by the minimum of two natural 
scales: (i) the localization length of non interacting parti- 
cles at zero energy C* = {h^ /VQam'^y^^ ^^'^^ and (ii) The 
healing length of the condensate = ^/\/m(/i + Vq). 
The case of smooth disorder, where a is the largest scale 
is conceptually somewhat simpler, and we shall therefore 
start the analysis from this regime. Later we will show 
that the rough disorder limit can be treated in an anal- 
ogous way. 



A. Smooth disorder 

In a smooth potential, the condensate has apprecia- 
ble amplitude only where the potential dips below the 
chemical potential. These regions define the superfluid 
puddles as illustrated in Fig. 1(a). Josephson coupling 
between puddles is induced by tunneling under the po- 
tential barriers separating them. 

Using Eq. (3) and the WKB approximation (see ap- 
pendix B) we obtain the coupling constants: 



J 



(4) 



Here Xi denote the two edges of the barrier, defined by 
V{xi)— jii = and = {dV/dx)^^^,. Hence the Joseph- 
son coupling is composed as a produ ct of three random 
variables J ~ ^1^2^, where yi = y^|F/| and T is the 
exponential factor. 

The macroscopic stiffness of the chain is determined 
by weak Josephson links that arise from atypically large 
barriers, much longer than the disorder correlation length 
a. For such barriers, the variables ^1, which depends 
on the left edge of the barrier, ^2 which depends on the 
right edge of it, and T which depends on the potential in 
the bulk of the barrier, are essentially independent. We 
can therefore obtain the distribution of each of the three 
variables separately in order to construct the distribution 
of J. 

To compute the statistics of the exponential factor in 
(4) for long barriers, we can split the integral to a sum 
on segments of size a on which the potential is approxi- 
mately constant 

1= dxy/2m{V{x) - /i) ^ cr ^ V^m{Vi - /i). (5) 
"^^1 i=i 

Here Vi > ji are independent random variables, dis- 
tributed as 



1 



(6) 
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= (27r) -^/^ J/T/Vo ^^^^ probability to find 

a potential V{x) > fi. At the same time, the proba- 
bility for a barrier of length / is p{l) = q^^^{l — Q/j,) ^ 
I Ing'^l exp(— I lng^|/), that is, the probability of having / 
consecutive segments with Vi > fi. Note that in the 2"^ 
(approximate) equality we moved to a continuous / while 
keeping the distribution normalized. 

The tail of the distribution of the sum / is dominated 
by large / (not by large V since the distribution of V 
has a much faster decay). We can therefore apply the 
central limit theorem to express / as a function of just two 
independent random variables I{l^r]) = a{Hil + rfSK^fl)^ 
where n and 5n are the average and standard deviation 
respectively of ^J2m{y — fi) over the distribution P^(V). 
77 is a gaussian variable with zero mean and unit variance. 

Given the distributions of / and 7^ it is a straight for- 
ward exercise to compute the probability distribution of 
the exponential factor T = exp {—I{r]^ 0/^)- limit 
of small T we find P{T) = AT^, with the exponent x 



and the pre-factor 
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as expressed in terms of the rescaled parameters n = 
V^2mVo/^(/i/Vb) and 5i<i = ^j2mVQ5'k{ii/VQ). A rigor- 
ous derivation of the distribution P{T) (see appendix C) 
gives essentially the same result as the more heuristic 
derivation presented above. 

The distribution of Josephson couplings may depend 
also on the distribution of the pre-factors yi = y^\V-\. 
Because the potential is a gaussian variable, i/i is dis- 
tributed as p{y) (X y dit small values of y. Consequently, 
the distribution of Josephson coupling at small values of 
J is given by J = ^1^2^ = (^/^o)(-^/^o)" with a = 
min{xi 1) • The power-law holds up to the cutoff scale 
determined by the pre-factor of (4). Using the typical 
value V ^ Vq/(7 we have l^o ^ ()A{h^ /rri)(yQ/cru). 

Note that some of the junctions in the effective model 
are actually formed with J > and therefore lie beyond 
the power-law distribution. However puddles separated 
by these strong junctions can be joined to make larger 
effective puddles in the same spirit of the real-space RG"^. 
If the charging energies are concentrated well below VLq^ 
which we shall see is the natural situation in experiments, 
then the process of removing strong links does not modify 
the exponent a of the distribution below l^o and we only 
obtain a normalized distribution for J G [0, VLq]: 



Po(J) = l^o'(c^ + l)(W)"- 



(9) 



When comparing to experiments it is important to take 
into consideration the physical length of the condensate 
and translate it to the number of junctions in the effec- 
tive Josephson array. By dividing the total length L to 



the average size of a junction and barrier we find the 
number = {L / (j)q^{l — q^). But as noted, some of 
the junctions lie above the cutoff VLq. After removing 
the strong junctions, as described above, we are left with 
Nq = dJP{J) = N^A {a + 1) junctions which fol- 

low the pure power-law distribution (9). This number 
turns out to be only slightly (< 10) lower than the total 
number of junctions for chemical potentials /i ^ Vq. 
In paractice it is therefore sufficient to take as the 
starting number of junctions. 

The charging energy of a puddle can be found using 
Thomas- Fermi approximation for the wavefunction of the 
puddle V^i(x) = ^/[jil^^V\x)y/u. The total number of 
particles is Ni = J^^'^ lipil'^dx and the total energy of the 
puddle is = I JJ^'* {i/jil'^dx where Xi^i and X2,i are the 
edges of the puddle. The charging energy is related to 
the chemical potential by /i = UiNi. Using the equations 
for Nj and Ej once can show that 



ON, 
u/Li 
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(10) 



which depends on the length of the puddle Li but not on 
the puddle's shape. 

We can compute the distribution of charging energies 
Ui in the array from that of the lengths Li of the super- 
fluid puddles. Since puddles are exactly complementary 
to barriers, the distribution of puddle sizes is the same 
as that computed above for the barrier lengths with q^ 
replaced by 1 — . The resulting distribution of charging 
energies is 



fo 



1^ r-^\ .-/o^o/t/+/o 



(11) 



where the cutoff scale Qq was estimated above. 

The distributions of the charging energies (11) and of 
the of the Josephson coupling constants (9) are precisely 
the stable solutions of the RG equations for the disor- 
der distributions^. The flow of these distributions, sub- 
stituted as initial conditions to the RG equations will 
therefore be greatly simplified: from the full functional 
flow to a flow of only two parameters a and /q. Such 
analysis will be described in section III. 

Finally we note that the results above are written in 
terms of the chemical potential, whereas the average den- 
sity is often easier to obtain from experiments. Given 
the gaussian distribution of the potential and using the 
Thomas Fermi approximation it is straight forward to 
obtain the relation between the two as 
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the same analysis, using Eq. (3), can be done to translate 
in-situ density profiles measured in experiments to a dis- 
tribution of Josephson coupling constants. The inset of 
Fig. 2 shows the exponent of the power-law as a function 
of the chemical potential derived from both the analytic 
and the numerical solutions. The approximate analytic 
result is seen to be in almost perfect agreement with the 
exact numerical calculation, though the former is rigor- 
ously controlled only for positive values of the chemical 
potential. 



III. RG FLOW AND RELATION TO 
EXPERIMENTS 



FIG. 2: Distributions of coupling Josephson couplings com- 
puted from exact numerical solution of the Gross-Pitaevskii 
equation for a condensate in a realistic disorder potential. In- 
set: the exponent of the power law fit of the low energy part 
of the distributions is compared to the analytic result using 
the WKB approximation. 



B. Rough disorder 

We turn to the case of a rough disorder potential, 
where a <^ C^^^h- Since the potential changes in space 
faster than the wave-function can respond, it is equiva- 
lent to a white noise potential, with the only important 
parameter being the combination D = VqCt. The dis- 
order strength D is directly related to the localization 
length at zero energy -'^^'-'^^ = {fi^ / Dm?y^^ and to the 
characteristic scale = /2mC,l. That is 1) = 

Since only the combination VqCf matters, while the in- 
dependent values of a and Vb are not important, the po- 
tential may be replaced by a smoothed effectiv e potential 
with cr^^ = min(C*,^^) and V^^ = Vo^/aja^ As be- 
fore we imagine filling this potential landscape, forming 
weakly coupled superfiuid puddles which grow with in- 
creasing chemical potential. In the puddles sets the 
smaller scale when ja > whereas is smaller for /i < 0. 
Under barriers on the other hand, C* is always the rele- 
vant length scale. Since we are interested in the vicinity 
of /i = where and are comparable we can apply the 
same analysis as outlined above for the smooth potential, 
taking cr ^ C* everywhere. This results in the distribu- 
tions (9) and (11) for the coupling constants. At high 
densities the exponent is a = 1 and it decreases together 
with the chemical potential and reaches the critical value 
a = when /i ^ h^/AmCl. 

In order to confirm the approximate analytic result we 
have solved the CP equation numerically for the case of a 
rough disorder potential. We find the ground state wave- 
function using imaginary time propagation and use it in 
Eq. (3) to obtain the distribution of Josephson coupling 
constants shown in Fig. 2(a) for a range of chemical po- 
tentials. The power-law behavior of the distribution at 
small values of J is seen clearly in the figure. Note that 



We have now derived an effective Josephson array 
model for bosons in a disorder potential that can be 
fed into the real-space RG framework of Ref. 4. The 
real-space RG consists of gradually eliminating sites with 
the largest charging energy Ui or Josephson coupling J^. 
The remaining sites are described by the same Hamilto- 
nian (2) with renormalized probability distribution of the 
charging energies and the Josephson coupling constants. 

The fiow of the effective disorder with decreasing en- 
ergy scale is described by a set of integro differential 
equations for the distributions of coupling constants. Re- 
markably the distributions of coupling constants Po{J) 
and Fo(/7), derived above from the microscopies, are pre- 
cisely self similar solutions to these equations'^. There- 
fore, when feeding these distributions as initial conditions 
to the RG equations the physics is fully determined by 
the fiow of the two parameters /o(r) and a(r) with the 
running RG scale 



rf/o(r) 
dr 

da{r) 
dT 



= /o(r) - /o(r) [a(r) + 1] 
= - [a(r) + 1] /o(r) 



(13) 



where T = log(f}o/^^)- 

The flow proceeds along the trajectories 

/o = a - In [a + 1] + e 



(14) 



shown in Fig. 3, where e labels the trajectory. The initial 
conditions /o(r = 0) and a{T = 0) of the fiow were 
derived above from the microscopic potential. The phase 
transition is crossed by changing the chemical potential 
or the disorder strength, which moves the initial point 
across the separatrix e = 0. 

In Fig. 3 we demonstrate the RG fiow derived for 
two realistic setups of ^^Rb condensates with disorder 
produced by a random speckle potential similar to Ref. 
1. Trap parameters are detailed in the figure caption and 
the disorder correlation length is taken to be the speckle 
size a ~ 0.26/im. The system is in the rough disorder 
regime in both cases. Experiments with atom-chip traps 
are in the opposite, smooth disorder, regime^^. But as 
we have shown, this leads to the same universal behavior 
on large scales. 
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FIG. 3: Real-space RG flow obtained from realistic setups 
with disorder produced by speckle potentials. The two sets 
of trajectories correspond to different values of the transverse 
frequency (jO± of the quasi one-dimensional trap potential and 
Id boson density p. uj± = 0.66 Khz{2 Khz) and p — 25(4) 
atoms per pm in the upper (lower) trajectories. The tran- 
sition (thick line) is crossed by varying the intensity of the 
speckle potential Vb as shown on the figure. The red boxes 
on the upper (red) trajectories mark where the flow should be 
terminated for condensates of total length 20, 100 and 500 pm. 



We also note that the specific structure of speckle 
potentials, non-gaussianity and spatial correlations, has 
known implications on Anderson localization. For ex- 
ample, emergence of pseudo mobility edges^^'^^. How- 
ever the exact structure of the speckle potential does not 
significantly affect the calculation of the Josephson el- 
ements. First, the Josephson elements are induced by 
tunneling of states under a long barrier and do not in- 
volve high energy states near the effective mobility edges. 
Second the exact probability distribution of V and in 
particular the fact that it is asymmetric does not enter 
the calculation in any important way. For our general 
framework to be applicable, the probability distribution 

of V should decay faster than e~>/^/^o at large positive 
V . If the decay is slower, then the weak links are not 
determined by rare long barriers but rather by the less 
rare high barriers. Using the exact structure of speckle 
potential in place of the Gaussian potentials used here, 
will result in a small change of the computed parame- 
ters a and /o, but not to a change in the form of the 
distributions. 

We can now also address the finite size of the sys- 
tem. For condensates of increasing length the RG flow 
should be terminated at decreasing energy scales, i.e. 
further along the flow, when all the elements in the 
effective Josephson array have been eliminated. The 
red squares on the flow trajectories in Fig. 3 demon- 
strate termination points corresponding to condensates 
of length 20, 100 and 500/im. The calculated values of 
physical observables, such as the superfluid stiffness or 
the compressibility^^, should be recorded at the termi- 



nation points and compared to the experimental mea- 
surements. To further characterize the critical point it 
would be interesting to study the coherence properties 
of the disordered condensate that can be extracted from 
interference experiments^^. This will require a general- 
ization of the theory of fringe statistics^^'^*^ to the case 
of disordered condensates. 



IV. PHASE DIAGRAM 

Above we have established a direct link between the 
microscopic model of bosons in a disordered potential and 
the renormalization group flow of the random Josephson 
array at strong disorder. This now allows us to com- 
pute the phase diagram in the space of microscopic pa- 
rameters, the dimensionless interaction and the disorder 
strength, at strong disorder. The transition line in that 
space is given by the set of models that map onto points 
on the separatrix as initial conditions for the RG flow in 
Fig. 3. 

The separatrix is given in terms of the RG flow pa- 
rameters by Eq. (14) with e = 0. We replace /o by its 
value as a function of the microscopic parameters (11) to 
obtain 

a{p) - In [a{jl) + 1] = -^u^ ln(l - g^). (15) 

where p = p/E^ and u = u/E^(^. The right hand side of 
the equation encodes the effect of quantum fluctuations, 
induced by the charging energies, which cause the bend- 
ing of the separatrix to positive values of a. By contrast, 
the mean field stiffness vanishes only for a < 0. 

The difference between the actual transition and the 
mean field approximation is most apparent when we plot 
the phase boundary in the space p versus u. In the 
regime of rough disorder, the disorder strength enters 
only through the energy E^. Therefore in this limit 
Eq. (15) charts a universal phase boundary in the space 
(/i, u). We can change p independently of the interaction 
u by tuning the density. However in the classical (Gross- 
Pitaevskii) solution the interaction enters only through 
the chemical potential, and so within this approximation 
the system becomes insulating below a critical chemical 
potential p^ ^ 0.47 independent of u. On the other hand 
the condition (15) gives the transition line 

p{u) = A* + O.S9u - OMu^ + O (u^) (16) 

showing a non trivial dependence on as a correction 
to the mean field result. The system is insulating for 
chemical potentials below the transition line. To obtain 
the formula (16) we have expanded both sides of Eq. (15) 
to quadratic order in p — p^. 

We stress that the quantum corrections stem from 
quantum phase slips, generated by the charging term in 
the effective hamiltonian (2). The RG flow accounts for 
such phase slips through decimation of sites with large 
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charging energies, which is accompanied by renormahza- 
tion of the local Josephson couplings. By contrast Bo- 
goliubov theory, being a quadratic expansion around the 
GP solution, cannot give rise to renormalization of the 
stiffness and therefore does not result in a correction to 
the mean field phase diagram^^. 

We can now find the phase diagram in the space of 
disorder strength versus interaction. The disorder in the 
rough limit is parameterized hy D = V^a = E^C,^. To 
obtain a universal phase diagram (i.e. independent of 
density and particle species) we turn to dimensionless 
parameters normalizing energies by the degeneracy scale 
Td = Ti?p^ /2m and length scales by the inverse density. 
Hence we use D = Dp/T^ and u = up/Td. 

As a first step we express Eq. (16) as a condition on 
the density rather than the chemical potential by using 
the equation of state (12) 

p ^ 0A9/U + 1.00 - 0.54u + O {u^) (17) 

where p = p(^. Finally, noting that u = u/p and D = 
1/p^ we can obtain the phase boundary in terms of u and 
D from (17) using a simple change of variables. This gives 
the asymptotic phase boundary for weak interactions 

\/^ ^ 1.70ii^/^ - 1.811^^/^ + O (li^/^) . (18) 

In order to obtain the phase boundary at larger inter- 
action strengths we take the exact form of the separatrix 
from Eq. (15) and use the exact equation of state (12) 
and solve the equation numerically. The phase bound- 
ary obtained in this way is shown in Fig. 4. It has the 
asymptotic form given by Eq. (18) at weak interactions. 
The mean field criterion a = 0, without the quantum cor- 
rections, gives only the first term in the expansion (18), 

i.e. V Dmf ^ 1.70'Lt^/^. This mean field transition line, 
previously found in Refs. [24] and [6] is also plotted on 
Fig. 4 for comparison with the actual transition. 

The phase boundary calculated using the mapping to a 
random Josephson array cannot be continued to arbitrar- 
ily strong interactions. Moving on the line given by Eq. 
(18) toward stronger interactions corresponds to moving 
the initial point of the RG flow on the separatrix shown 
in Fig. 3 from the point a = /o = toward higher values 
of a. This is tantamount to decreasing the disorder on 
the effective Josephson array, which makes the starting 
point of the RG flow gradually less controlled. Specifi- 
cally, at the dimensionless interaction strength > 1.5 
the exponent a reaches the value a ~ 1. This can serve 
as a characteristic value beyond which the distribution 
of Josephson links becomes narrow and where the real 
space RG approach is not justified a priori. 

To complement the phase diagram found in the weak 
interaction limit, we also obtain the phase boundary be- 
tween the superfluid and the insulator in the regime of 
weak disorder. This is done using the perturbative RG 




FIG. 4: \/D-u phase diagram. The diagram is given in terms 

of dimensionless disorder strength \/d — ^Dp/Td and the 
dimensionless interaction strength u = up/Ta^ where p is the 
density, u is the effective contact interaction, D is the measure 
of rough disorder defined in the text, and Td = p^ /2m is the 
temperature of quantum degeneracy. The phase boundary at 
weak interactions corresponds to the set of points mapped to 
the separatrix of the RG flow. The dashed line above the true 
transition marks the points where the classical theory predicts 
a transition^"^ based on vanishing of the mean field stiffness. 
The phase boundary at strong interactions is obtained from 
the weak disorder theory of Giamarchi and Schultz^(GS). 

flow of Giamarchi and Schulz^ 



2| = (3-2A-')P 

Here we used the convention of Ref. [9] for the Lut- 
tinger parameter K. By the dividing the two equations, 
the separatrix of the flow is easily found to be given by 
b = 6{K-Ko)-4.ln{K/Ko), where Kq = 2/3 is the value 
of the Luttinger parameter at the fixed point. Now, us- 
ing the relation of the bare Luttinger parameter (i.e. K 
of the clean system) to the microscopic interaction in the 
relevant regime ^ 1 + 4/'u,^^ we obtain the approx- 
imate phase boundary ^ {8 — u)/\2\f2. Of course 
this approach is justified only as long as the disorder 
strength D is smaller than 1. The global phase diagram 
inferred from both the strong and weak disorder limits is 
plotted in Fig. 4. 

V. CONCLUSIONS 

In this paper we derived an effective quantum Joseph- 
son array model starting from a realistic microscopic 
model of bosons in a one dimensional random potential. 
The distributions (9) and (11) of the effective coupling 
constants obtained in this way are precisely the stable so- 
lutions of the real-space RG found in Ref. 4. Their flow 
is determined by a strong randomness fixed point, which 
controls a quantum phase transition between a superfluid 
and insulating phase. 
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The ab-initio mapping from microscopic models of 
bosons in a random potential to precise initial conditions 
for the RG flow allows to make quantitative predictions 
and thereby can facilitate experimental detection of the 
new critical point. This can be done for example by mea- 
suring the finite size scaling behavior of the phase corre- 
lations in interference experiments. The results can be 
directly checked for consistency with the finite size scal- 
ing implied by the RG flow near the critical point. 

Bridging the gap between the microscopic physics and 
the RG flow of the Josephson array model allowed us 
to predict a phase diagram (Fig. 4) in the space of mi- 
croscopic parameters, disorder and interaction strength, 
which is valid in the regime of weak interactions. The 
transition line is modifled compared to that previously in- 
ferred from mean fleld (Gross-Pitaevskii) theory^^'^^ and 
from estimates based on typical values of the Josephson 
coupling compared to interaction strength^. 

Note that the critical point in this regime is different in 
nature from the weak disorder transition considered by 
Giamarchi and Schulz^, shown as a separate transition 
line in the strong interaction range of Fig. 4. At the 
strong randomness flxed point rare events in the form of 
weak links that effectively cut the chain play a central 
role, while these are completely neglected in the weak 
disorder theory. As shown in Ref. [11], such events lead 
to a transition at a non universal value of the phase- 
correlation decay exponent that is always smaller than 
the universal value of the exponent predicted by the weak 
disorder theory. It remains an interesting open question 
how the transition interpolates between the two limits. 
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Appendix A: Superfluid Stiffness 

The Gross-Pita evskii energy functional of the wave- 
function tlj{x) = ^/ p{x)e^'^^^^ can be written as: 



Jo 







dx(^-^^^V{x)p^lup'' 
Amp 2 ^ 



(Al) 



where A is a Lagrange multiplier that can be used used 
to enforce a phase twist ^ = dx^f along the length of 
the condensate . For <l> = the ground state has uniform 
phase and can be chosen to be real i/jq = ^/pO' The energy 
of this state is cq = Eo[po]. 

The superfluid stiffness is proportional to the quadratic 
change of the energy with a small twist $ 



L 



(A2) 



To compute ps we write the ground st ate wave-fun ction 
in presence of the twist as V^<^>(x) = y^po + 8p{x)e^'^^^\ 
and expand the energy to quadratic order in the changes 
5p and dx^ 



E^=Eo^O{Sp^)- 



Xdx^) . (A3) 



We now note that the phase twist ^ led to a propor- 
tional phase gradient dx^, while the change in the local 
density must be quadratic Sp cx This is because p 
is even under time-reversal whereas ^ is odd. Therefore 
the change in energy due to distortion of the density by 
the twist is proportional to and does not contribute 
to the stiffness. 

Minimizing (A3) with respect to the phase gradient we 
have 



dx(p{x) 



m A 



(A4) 



and A = 



is obtained by imposing the con- 



straint. By substituting back in Eq. (A3) we flnally ob- 
tain the superfluid stiffness 



L- 



1 



(A5) 



Appendix B: Josephson Coupling 

Using the above equation for the SF-stiffness, we may 
calculate the effective Josephson coupling of two neigh- 
boring SF puddles. The energy of the coupled puddles is 
E ^ Jcos{^i — $2)- Therefore, it can be related to the 
SF stiffness of the system in the region between the two 
puddles by 



J: 



12 



(Bl) 



where xi {X2) is the left (right) edge of the barrier be- 
tween the two SF puddles. We are left with the problem 
of flnding the wave-function in the region between the 
two puddles. This is done in a way similar to the calcu- 
lation in Ref. 27. 

since the wave-function amplitude is small under the 
barrier, the interaction term is negligible. The 

wavefunction follows a linear Schrodinger equation and 
it can be approximated using the WKB approximation 



^o(^) 



^1 ^- dyKiy) 



Co 



^/K{x) 



^/K{x) 



<3E>=0 



with n{x) = y^[p — y{x)] and p the chemical poten- 
tial. The effect of the puddles enters through the pre- 
factors Ci and C2. These are determined by the match- 
ing conditions of the WKB wave-function with the wave- 
function inside the puddles. 
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In the regime where the heahng length of the puddles 
is smaller than the size of the puddles, we may use 
the Thomas-Fermi (TF) wave-function inside the pud- 
dles V^o(^) = \f^^^^- However, the TF-approximation 
breaks down at distances closer than ~ ^/j, from the edge 
of the puddles. In this narrow region the potential can 
be taken to be linear, and since the amplitude is small 
we can again drop the interaction term. The solutions 
to the Schrodinger equation in this region are Airy func- 
tions, which we can match with the WKB wavefunction 
and with the TF wave-functions on the two sides. This 
gives the constants: 



a; 



2u Ai [2-2/3] 



0.4 



\V'{x.)\ 



. (B3) 



Having obtained the wave-function V^o(^) we are in po- 
sition to compute the superfluid stiffness using (Bl). Be- 
cause i/jQ is exponentially suppressed in the middle of the 
barrier the integrand IV^ol"^ is strongly peaked suggest- 
ing the use of a saddle point approximation to evaluate 
the integral. We write the integral as / dx\ilJo{x)\~'^ = 

j dxK{x)ef^''^ ^ /(ixA>:(xo)e'^^^o)+i'^"(^o)(^-^o)', with 
the obvious definition of f{x). We approximate it using 
a saddle point approximation around xq which is the po- 
sition of the maximum of f{x). By differentiating f{x) 
one can show that xq satisfies f"{xQ) = —2hz'^{xo) and 



e 2 



. Using these in the 



saddle point approximation, the Josephson coupling is 
found by gaussian integration 

^ ^'ciC2e-^-''^"'^^"\ (B4) 



Jl2 = 



/TT m 



J is the product of the tunneling coefficient, denoted by 
T, and a non-trivial pre-facor. 



Appendix C: Distribution of the tunneling coefficient 

The tunneling coefficient T is the exponential factor in 
eq. (B4) 



T = e 



(CI) 



We shall approximate the integral in the exponent as 
a sum on segments of size a, the correlation length of 
the potential, and take V{x) = F^, a constant on each 
segment 



rX2 

J X\ 



1= I dxJ^{V{x)-f,) 



(C2) 



i=l i=l 

The distribution of the potential in each segment is given 
by the conditional probability 



P^{V) = PiV\V > ^) = 



1 



(C3) 



where = P(y > ji) is the integrated probability for 
the potential to be higher than the chemical potential fi. 
The number of summands / = {x2 — Xi)/cf is a random 
variable in itself with the distribution 



p{i) 



q'-\l - q^) « \\nq^\ exp(-|lnqJO (C4) 



In the last equality we have taken the continuum limit of 
I. This is not expected to affect the distribution of T in 
the limit of small T, since this tail is controlled by large 
values of /. 

We can now derive the probability distribution of the 
sum /. Formally it is given by 

P{I) = / dlp{I\l)p{l) (C5) 
Jo 

where p{I\l) is the conditional probability of / given a 
barrier length 1. Since / for a given / is a sum of in- 
dependent and identically distributed random variables 
the characteristic function is simply related to the 
characteristic function of each term Zi as ^i{t) = (l)z{ty. 
The conditional distribution p{I\l) is given by the inverse 
fourier transform of the characteristic function 



Pirn 



dtexp{-itl)4>ziity 



(C6) 



Putting together Eqs. (C4), (C5) and (C6) we get 



p{I) 



2 roo poo 

— / dl dte-'''Mty\^M^\e~ 

^TT Jo J-oo 



27r 
IW^I / 

27ri 



dte 



-iti 



\\nq^\-\n{4>,{t)) 



ZOO 



I 

J —1 



1 



llng^l -ln(0^(-zr))* 



In the last line we changed variable r = it, and integrated 
over the imaginary axis in the complex plane. 

We can close the integration contour at infinity over 
the right semi-circle, where the integrand decays expo- 
nentially, without changing the integral (see Fig. 5). The 
contributions to the integral then come from the poles at 
the points having a positive real part 



p{i) = J2 ^~ 



^^Res 



\nq^ + ln [(l)z{-iT)] 



(C7) 

The points are zeros of the denominator and therefore 
obey the equation (j)z{—iTi) = 

We shall now prove a theorem regarding the solutions 
of the above equation. This theorem will enable us to 
predict the tail of the distribution of T. 

theorem 1. The pole condition (j)z{—ir) = 1/q^ has a 
real solution tq > 0. All other solutions Ti satisfy Reri > 
To- 
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Proof. First, we prove the first part of the theorem. Us- 
ing the definition of the characteristic functions with real 
argument r we get the pole condition 



f 

Jo 



(C8) 



Note that the lower limit of the integration is zero and 
not — oo since p{z < 0) = 0. is a probability, hence the 
R.H.S. of (C8) is larger than 1. The integral on the L.H.S. 
converge for any real r since p{z) decays faster then any 
exponential. By differentiation one can show that this 
integral is a monotonically increasing smooth function of 
r. For r = it is 1. Therefore, there is always some real 
To > which satisfies the pole condition. 

We prove the second part of the theorem by contra- 
diction. Let us assume that there is a complex pole at 
Ti = a -\- ib with a < tq. Both and tq obey the pole 
condition, and so (pzi—iTo) = = (j)z{—iTi). Using 

the definition of the characteristic functions this equality 
can be written as: 



0= / d^(e^°^ -e^^^)p{z) 



dze""' (e^^o-a)z _ ^ibz^ ^^^^ 



(C9) 



Now, by our assumption that a < r the real part of the 
expressions in the brackets above is positive (it may be 
zero on a discrete set of points if a = tq). Since all other 
factors in the integrand are non-negative, we conclude 
that the integral cannot vanish, which contradicts our 
assumption. Therefore we must have a = Reri > tq. □ 

The distribution of T as T ^ is determined by the 
tail of the distribution p{I) and it is therefore dominated 
by the pole tq with smallest real part. This pole must 
be real according to the above theorem. It follows, by 
change of variables, that the distribution of T behaves 
as a power-law in the limit of small T. That is, P{T) = 
j^j'ro-i^ where A is the residue of the pole at tq. 

To obtain the numerical value of the exponent we have 
to find a real solution to the pole condition. This is done 
by expanding the log of the characteristic function in 
powers of r 



ln[(j)i{T)] = 22 



(CIO) 



where /^^ is the n*^ cumulant of the distribution p{zi). 
The first two cumulants are hzi = (zi) = Zi and = 
(zf) — (zi)'^ = 5zf. In general, the expansion has some 
finite radius of convergence. We will come back to this 
point later. 

Plugging the expansion in the pole condition (j){—ir) = 
gives the infinite series equation 



(Cll) 



I mi 



Rei 



FIG. 5: The contour integration. The contour integration in 
the complex plane over the right semi-circle. The integral 
collects all the poles in the right half plane 
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FIG. 6: Zeros of the truncated polynomial. The zeros of the 
truncated polynomial (31^* order) in eq. (Cll) which corre- 
spond to poles in the contour integration, as calculated for 
/i = 0. Spurious zeros cluster along a circle, which is the ra- 
dius of convergence of the cumulant expansion. The real zero 
is inside the circle. 



To approximate the real zero we truncate the series 
and find the zeros of the resulting polynomial. The 
zeros of the truncated polynomials have the following 
properties'^: 

1. If the zero is within the radius of convergence, the 
truncated polynomial will converge to the true zero 
as the length of the polynomial is increased. 

2. If the zero is outside the radius of convergence, it 
will not be found. 

3. Zeros of the polynomial which are not really zeros 
of the infinite series will cluster along the radius of 
convergence. We call these zeros spurious zeros. 

Since the real zero that we look for is generally within 
the radius of convergence, we may truncate the series to 
approximate its value. We keep only the first two terms 
of the cumulant expansion, which gives us a 2"^ order 
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polynomial, with the positive solution 

To , (Clij 

and the residue, which is also the pre-factor of the power- 
law tail 

(|lng,|-ln(0,(r))),^ " "^^'^T^^IIifS^^ 

(C13) 



Taking higher orders of the polynomial changes the po- 
sition of the real pole only slightly. In fact tq is approx- 
imated by (C12) with less than 0.03% error. All other 
zeros of the polynomial cluster along the radius of conver- 
gence, as expected. This is demonstrated for expansion 
up to 31^* order in fig. 6 for /i = 0. 
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